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Abstract. We propose that an onset of glass transition can be defined as the 
point at which a supercooled liquid acquires the stress relaxation mechanism of a 
solid glass. We translate this condition into the rate equation for local relaxation 
events. This equation simultaneously gives two main signatures of glass transition, 
stretched-exponential relaxation and the Vogel-Fulcher law. The proposed theory 
quantifies system fragility in terms of the number of retarded local relaxation 
events and reproduces its correlation with the non-exponentiality of relaxation 
and bonding type. 
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If a liquid is cooled down fast enough, it forms glass. At the onset of glass 
transformation range, a liquid qualitatively changes its properties, and the two 
main features that distinguish it from a high-temperature liquid are the stretched- 
exponential relaxation (SER) and the Vogel-Fulcher (VF) law. When a perturbation, 
in the form of stress or external field, is applied to a liquid near glass transition, a 
relaxing quantity q(t) decays following a universal SER: 



where < j3 < 1. This behaviour is seen in many systems, and is considered a 
signature of the "glassy" relaxation ^ El ID IS] • Another universal feature of this 
regime is that viscosity, or relaxation time r, follows non-Arrhenius dependence, which 
in most cases is well approximated by the VF law: 



where A and To are constants [5]. 

A substantial amount of research in the area has revolved around the origin of 
these two anomalous, yet universal, relaxation laws. A successful theory of the glass 
transition, as widely perceived, should provide a common justification for Eq. (1) 
and (2) 0Ej- Recently, the need for a theory of the glass transition to give Eq. 
(1) and (2) simultaneously has been reiterated on the basis of the close relationship 
between (3 and r: it has been found that (3 is invariant to different combinations of 
pressure and temperature that hold r constant [Jj. It has therefore been suggested 
that this correlation should constrain any theory of the glass transition, in that if a 
given formalism gives Eq. (2), it should also be able to give Eq. (1) 0. 

Several decades ago, Goldstein proposed (S] that at glass transition, flow becomes 
dominated by potential barriers which are high compared to thermal energies, whereas 
at high temperature, the opposite is true, barriers are much smaller than thermal 
energies. Hence it can be argued that while in a liquid at high temperature, local 



q(t) K exp(-(t/r)' 3 ) 
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stress is relaxed on the timescale of microscopic trajectory reversal times, in the 
supercooled regime, local regions can support a finite stress (i. e. maintain local 
structure unchanged) on timescales that are considerably larger. This has opened 
the possibility to discuss the stress relaxation mechanism in a liquid approaching 
glass transition. However practical realizations of this approach, in particular the 
relationship between the stress relaxation mechanism and Eq. (1) and (2) have 
remained elusive. 

Following this approach, we consider that as temperature decreases, liquid 
acquires a "solid-like" ability to support local stresses on timescales that considerably 
exceed trajectory reversal times. Our main proposal is that a liquid near glass 
transition also starts to redistribute local stresses in a solid-like manner. Hence the 
onset of the glass transition can be robustly defined as the point at which a liquid and 
a solid glass under stress begin to redistribute local stresses in the same way. In other 
words, we propose that the onset of glass transition is the point at which the liquid 
acquires the stress relaxation mechanism of solid glass. We show that this condition is 
sufficient to recover both anomalous relaxation laws, Eq. (1) and (2). We also show 
that the proposed theory gives a simple definition of a system's fragility in terms of the 
number of local relaxation events induced by external perturbation, recovers fragility 
plots and predicts correlations of fragility with (3 and the nature of the chemical bond. 

How do local stresses redistribute in glass under pressure? In the same paper |S] , 
Goldstein considered this question: he argued that because a local region supports 
less stress after the relaxation event than before, all other local regions support more 
of the external stress after the event than before. Generally, increasing stress on 
other local regions makes their relaxation more difficult. The increased stress they 
need to support is aligned along the direction of external pressure. On the other 
hand, local relaxation paths for local events with the smallest barriers are generally 
oriented at random relative to the external pressure, because they are defined by the 
symmetry of local ordering 9 (in Goldstein's terminology, local reaction paths are 
"non-concordant" to the external stress ;8 ). As a result, activation barriers increase 
for later local relaxation events. 

In what follows, we consider that relaxation proceeds by local relaxation events 
(LRE). In the literature, these jump or flow events have been given different names (for 
review, see Ref. JO]); m this discussion we borrow the term LRE from our previous 
studies of relaxation in glasses [HI El El EH- In glass under high enough pressure, a 
LRE involves several localized atomic jumps which include breaking old bonds, forming 
new ones and the subsequent relaxation of the local structure [5]. An animation of a 
LRE in Si02 glass is available in the electronic form of Ref. [5]. Each LRE carries a 
microscopic change of a macroscopic relaxing quantity, e.g., volume. By considering 
the dynamics of LRE and their coupling to structural rigidity of glass, it has been 
possible to explain several interesting aspects of glass relaxation, including the origin 
of slow relaxation [111 El an d the origin of temperature-induced densification in the 
pressure window, centered at the rigidity percolation point J3| • 

We introduce LRE as local relaxation "quanta" which a liquid uses to adjust to 
external perturbations. Each LRE carries a microscopic change of a liquid's relaxing 
quantity (i. e. volume, external stress etc). In a high-temperature liquid, a LRE 
is an atomic jump from the surrounding "cage" , followed by local relaxation. As 
temperature decreases, atomic rearrangements become more cooperative, due to the 
need to cross higher activation barriers. In this regime, each LRE is associated with 
the transition over the activation barrier in the Goldstein picture of activated flow [BJ • 
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In the supercooled regime and below, LRE, induced in different parts of a system, have 
different relaxation times, i.e. they are dynamically heterogeneous |14j . as discussed 
below in more detail. 

We now derive the rate equation for LRE in a liquid at the onset of glass transition. 
For this, we use our main proposal that at glass transition, liquid acquires the stress 
relaxation mechanism of solid glass. So first we find how to express the stress relaxation 
mechanism in glass in mathematical terms. In particular, we find how activation 
barriers for LRE change as a result of redistribution of local stresses. Let N be the 
total number of relaxing units in the structure. Under external (hydrostatic or shear) 
stress P, each unit supports stress po such that P = poN. Since, as discussed above, 
after relaxation, a local unit supports stress p\ < po, the stress on other local regions 
is Pi = (P — Pin c )/(N — n c ), where n c is the current (instant) number of LRE induced 
by external perturbation. If n = n c /N <c 1, p% — Pa + {po — Pi)n. It has been 
argued that the main contributor to the activation barrier V is elastic energy [TU] , 
Hence the increase of V is proportional to the increase of work needed to overcome 
the barrier created by elastic force due to additional stress Ap = p2 — pa = (pa —pi)n. 
So V oc Ap oc n for small n: 

V(n) =V + Vm (3) 

where Vq is the energy barrier in an unperturbed system, and V\ is defined such that 
V(n r ) is the maximal energy barrier, where n r is the total number of LRE caused by 
an external perturbation, n(t) — > n r as t — ► oo . 

Note that Eq. © can not be applied to a liquid above glass transition, because 
at high temperature externally-induced stresses are quickly removed by thermal 
fluctuations, and redistribution of stresses between different local regions does not 
take place. More precisely, at high temperature, stress relaxation mechanism, as 
described by Eq. J2J, only exists on short microscopic timescales. On experimental 
timescales, over which Eqs.lQJ and are measured, V is independent of n at high 
temperature. 

We now apply Eq. Jjyl to the liquid approaching glass transition. First, the rate 
of LRE, 4j, depends on the event probability, exp(— V/kT). According to our main 
proposal, the onset of the glass transition is defined as the point at which a liquid 
acquires the stress relaxation mechanism of a solid glass; hence V is given by Eq. 
(0. Second, because an external perturbation induces a finite number of relaxation 
events, n r , the rate of LRE should also have a saturation term to reflect the depletion 
of LRE. The most natural choice for the saturation term is linear —an dependence, 
which reflects the fact that relaxed events are removed from further dynamics. This 
is analogous to, for example, the process of nuclear decay, in which the decay rate 
decreases linearly with the number of decayed nuclei, An/ At oc — n. Hence, using Eq. 
J3J and assuming that in a liquid Vq <C kT, we write 



— = exp(— Cn) — an 

where C = Vi/kT and t is re-scaled &s t t/to, where to is the characteristic 
relaxation time, a is defined from the condition that = when n = n r , giving 
a — exp(— Cn r )/n r : 

= exp(— Cn) — — exp(— Cn r ) (4) 
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Before solving Eq. we note that the only assumption in its derivation is 

that the increase of the activation barrier is linear, Eq. @. We argued that this is 
justified for n <C 1, and here we note that Eq. is consistent with experiments in 
common glasses. If the saturation effects are small and the second term in the right 
part of Eq. is ignored, n oc ln(t + 1 ). Thus linear expansion gives logarithmic 
dependence of n, which is consistent with the logarithmic relaxation of macroscopic 
properties (e.g., volume) of SiC>2 and GeC>2 glasses under pressure jTHj. In addition 
to the qualitative agreement, Eq. J3J also gives a quantitative agreement with the 
experiment, as is found by the calculation of the slope of logarithmic relaxation of 
volume under pressure 

Eq. Q has two parameters, C and n r . When Cn r <C 1 (and hence Cn <C 1 since 
n < n r ), the right part of Eq. (@J becomes 1 — n/n r , leading to the usual Debye- type 
relaxation. This takes place when either the temperature T oc 1/C is high, or when 
n r is small. At Cn r w 1, one expects the onset of non-exponential relaxation. This 
sets the scale for the non-exponentiality temperature T n : 

kT n w V\n r (5) 

Below we show that in the T « T n regime, Eq.(4) gives SER and the VF law, Eqs. 
(1) and (2). Note that T n is higher than the glass transition temperature T g . T g is often 
defined from the condition of viscosity reaching some large value, corresponding to the 
relaxation time exceeding the time of experiment, whereas T n defines the preceding 
temperature regime at which relaxation becomes non-exponential and non-Arrhcnius. 

First, we solve Eq. Q for a wide range of parameters (C, n r ) that satisfy condition 
Cn r > 1, i.e. define the non-exponential regime of solution of Eq. J3J. Remembering 
that n(t) — ► n r as t — > oo, we fit the solution to n = n r (l — cxp(— t/r)"). We note that 
this form of SER and Eq. (@J contain two parameters each, which suggests that if a 
good fit exists, it is not accidental, but probably reflects the involved physics. Figure 
1 shows that fits of the solution to SER are very good. We find that this is the case 
in the wide range of parameters (C, n r ), except when Cn r ^$>l. 

An important observation from Figure 1 (see the legend) is that (3 decreases as 
C'n r increases. Hence non-exponentiality can increase as a result of either increase of 
n r , or decrease of temperature T oc 1/C. This is consistent with many experiments 
in which (3 decreases with T [3 El El- We will return to this point below. 

Second, for several different n r , we solve Eq.(4) as a function of C, and fit the 
solution to the SER form above to obtain r. We plot the solution as a function of 
l/(n r C)= T/T n , where T n is defined from Eq. (0). We find that the dependence of 
relaxation time r on T/T n collapses on the curve ]n(r/n r )=f(T/T n ). We also find 
that f(x) can not be represented by the Arrhenius-type dependence oc 1/x, however 
a good fit is obtained if 

In — = fl i = «i ^ 

n r l/(Cn r )—CL2 T/T n — a,2 

where a\ and are constants. This is the form of the VF law, Eq. |J2J. We find 
that a good fit to Eq. exists in both T > T n and T < T n regimes, with ln(r/n r ) 
spanning over 15 decades (see Figure 2). 

The transition from Eq. © to Arrhenius dependence directly follows from the 
condition that solution of Eq. Q is exponential, Cn r <C 1. Since 02 < 1 (see the 
legend to Figure 2), this condition means l/(Cn r ) 3> 02, and dependence becomes 
Arrhenius. If, on the other hand, Cn r <C 1 does not hold, relaxation is VF-type. 
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Figure 1. Solid lines are the solutions of Eq. (4) for several pairs of parameters 
(n r ,C): 1 - (0.001,1000), 2- (0.001,4000), 3- (0.1,10). 4 - (0.1,50), 5 - (1,3). 
Dashed lines are the least-square fits to SER, giving the following parameters of 
t): 1- (0.93,0.0011), 2 - (0.63,0.0031), 3- (0.94,0.114), 4 - (0.55,0.478), 5 - 
(0.71,2.05), respectively. For each value of n r , the solution of Eq.(4) for n has 
been divided by n r so that 0< n/n r <1. 
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Figure 2. The solid line is the solution of Eq. (4), fitted to SER to obtain r. 
The dashed line is the fit to Eq. (6), with a\ =0.141 and a2=0.116. The insert 
shows the fit in the T < T n regime; ai =0.288 and a2=0.017. 

We now suggest that the proposed theory clarifies the origin of system "fragility" , 
its correlation with the non-exponentiality of relaxation E] and the nature of the 
chemical bond. 

First, n r , a parameter in Eq. quantifies the overall atomic motion in a system 
due to external perturbation that comes in addition to thermally-induced motion. At a 
given temperature, n r depends on the magnitude of external perturbation and, more 
importantly on the system's ability to resist structural changes at the microscopic 
level. This ability has been termed a system "fragility", and constitutes the basis 
of fragility plots, which are essentially plots of Eq. JSJ with a varying parameter 
that measures the deviation from Arrhenius dependence; the larger this deviation, the 
larger the fragility 
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Figure 3. Decrease of j3 with n r (fragility), at various values of C. At constant 
n r , (3 increases with T oc 1/C. 

Qualitatively, a "strong" system has a built-in resistance to temperature-induced 
structural changes, whereas the structure of a "fragile" system is easy to disrupt 5 . 
In our picture, this means that a strong system responds to external perturbation 
with little retardation and smaller n r , whereas relaxation in a more fragile system 
involves a larger number of LRE, required to come to equilibrium with new external 
conditions. In other words, a more fragile system is more retarded in terms of larger 
n r . 

Hence our picture offers the quantification of fragility in terms of n ri and we 
can immediately recover fragility plots in this approach. As discussed above, small 
n r corresponds to a stronger system. At a given temperature, small n r results in 
l/(Cn r ) 3> a 2 , and Eq. © becomes Arrhenius. As n r increases, corresponding to 
a more fragile system in our picture, l/(Cn r ) S> a 2 does not hold, and relaxation 
becomes progressively non- Arrhenius (see Eq. JB}). Another way of discussing this 
effect is to note that n r is proportional to T n (see Eq. |5j); hence larger n r in a fragile 
system corresponds to higher T n . Physically, this means that the increase of the degree 
of system's retardation, quantified by n r , requires a higher temperature to remove this 
retardation by a more efficient equilibration and make the system relax exponentially. 

Second, we find that the proposed picture reproduces the relationship between 
fragility and non-exponentiality. Experimental data of more than 70 systems 4. show 
that [3 decreases linearly with fragility. In our picture, fragility is defined by n r , and 
in Figure 3, we plot (3 as a function of n r for different values of C. It is indeed seen 
that decreases with n r , reproducing the experimental correlation well. Note that at 
a given n r , higher T oc 1/C results in the increase of j3 (see Figure 3), in agreement 
with experimental observations El El- 

Finally, one can discuss how the proposed picture relates a system's fragility to its 
microscopic parameters. Other conditions being equal, one expects that an external 
perturbation induces generally larger n r in a system with ionic bonding as compared 
with a system with covalent bonding. In the covalent case, an atomic pair lowers its 
energy through sharing electrons between two atoms, resulting in a binding energy 
as high as several eV, and a LRE necessarily involves breaking these stable electronic 
configurations (breaking "covalent bonds") with associated high energy cost. In the 
ionic case, atomic rearrangements can proceed without a change in the electronic 
state of the atoms. As a result, activation barriers generally increase with covalency 
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of bonding. Since, as discussed above, fragility increases with n r , one readily predicts 
that covalent systems should be generally stronger and ionic systems should be more 
fragile, in good agreement with experimental results [SJ- Other factors, in addition to 
the nature of the chemical bond, may also affect n r , including for example, the ratio 
of ionic radii. 

Before concluding, we make three remarks. First, it is important to note that 
Eq. Q yields SER and the VF law simultaneously, suggesting that LRE dynamics 
are behind both anomalous "glassy" relaxation laws that kick in at the onset of the 
glass transition. This clarifies an open question of why the relaxation function is non- 
exponential at temperatures at which the relaxation time is non-Arrhenius [Hj : in our 
theory, larger n r increases the non-exponentiality of relaxation (see Figures 1 and 3) 
and, at the same time, increases departure from the Arrhcnius relaxation as follows 
from Eq. ©. 

A second related point is that the relationship between SER and the VF law 
has recently been reiterated: it has been discovered that j3 is invariant to different 
combinations of pressure and temperature that hold r constant |7| . It has therefore 
been suggested that this correlation should constrain any theory of the glass transition 
jZj. In our theory, temperature and pressure define parameters C and n r in Eq. J1J). 
C and n r , in turn, unambiguously define (3 and r. Due to the monotonous character 
of the solution of Eq. we find that only one value of (3 corresponds to a given r. 
In other words, we find that j3 is invariant to different combinations of pressure and 
temperature that hold r constant, satisfying the experimental result [Jj. 

Finally, we note that in order to derive Eq. (1) and (2), we considered a system 
under external stress, which allowed us to discuss the stress relaxation mechanism 
near glass transition, Eq. II) . At the same time, Eq. (1) and (2) are observed 
in supercooled liquids in the absence of pressure as well, from decay of correlation 
functions. This behaviour directly follows from the considered situation of the system 
under stress, by applying the fluctuation-dissipation theorem. 

In summary, we proposed a new simple way of defining the onset of glass 
transition: a liquid enters glass transformation range when it begins to redistribute 
local stresses in the same manner as solid glass. We showed how this condition 
simultaneously gives two main signatures of "glassy" relaxation, the stretchcd- 
exponential relaxation and the Vogel-Fulcher law. Consistent with recent experiments, 
we found that in the proposed theory, [3 is invariant to different combinations 
of pressure and temperature that hold t constant. We have discussed that the 
proposed theory offers the definition of system's fragility in terms of the number of 
local relaxation events, and recovers experimental correlations of fragility with non- 
exponentiality of relaxation and the nature of the chemical bond. 
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